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ABSTRACT 

Aims. Chemical databases such as the UMIST Database for Astrochemistry (UDFA) are indispensable in the numerical modeling of 
astrochemical networks. Several of the listed reactions in the UDFA have properties that are problematic in numerical computations: 
Some are parametrized in a way that leads to extremely divergent behavior for low kinetic temperatures. Other reactions possess 
multiple entries that are each valid in a different temperature regime, but have no smooth transition when switching from one to an- 
other. Numerically, this introduces many difficulties. We present corrected parametrizations for these sets of reactions in the UDFA06 
database. 

Methods. From the tabulated parametrization in UDFA, we created artificial data points and used a Levenberg-Marquardt algorithm 
to find a set of improved fit parameters without divergent behavior for low temperatures. For reactions with multiple entries in the 
database that each possess a different temperature regime, we present one joint parametrization that is designed to be valid over the 
whole cumulative temperature range of all individual reactions. 

Results. We show that it is possible to parametrize numerically problematic reactions from UDFA in a form that avoids low temper- 
ature divergence. Additionally, we demonstrate that it is possible to give a collective parametrization for reaction rate coefficients of 
reactions with multiple entries in UDFA. We present these new fitted values in tabulated form. 
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1. Introduction 

Astrochemistry is a field of increasing scientific attention. Over 
the last three decades, a lot of effort went into understanding 
how interstellar chemical species are formed, destroyed, and 
how interstellar chemistry works in detail. F or a review of re - 
cent as t rochemical progress see for insta n ce Lis et al.l d2005); 
iHerbsd (l2006l) : iLemaire & Combesl (I2007I) : Ivan der Takl (E008). 
To model interstellar chemistry, reliable laboratory data on the 
respective reaction rate coefficients are of utmost importance. 
A few comprehensive databases of astrochemical rea ctions are 
available today, such as UDEAQ (Wood all et"ai1l2007l) . the Ohio 
database OSL0 (e.g. jGarrod et al.ll2008l) . and NIST Chemistry 
Webboolfl There are also efforts to pool all available reac- 
tion data into a uni fied database ( KIDA: Kinetic Database for 
Astrochemistry j3 (Wakelam 2009). Results from astrochemical 
laboratory experiments are collected in these databases in an 
indispensable effort to make these results known to the wider 
community. They represent the interface between chemistry and 
astrophysical modeling. In the following we will mainly con- 
centrate on UDFA06, which contains parametrized reaction rate 
coefficients for more than 4500 reactions, involving 13 differ- 
ent elements and 420 different chemical species - far more than 
have actually been observed in interstellar space. 
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UDFA (and most other databases as well) tabulate two-body 
reactions together with a set of fit parameters a, (3, 7 in the form 
k(T) = a(T/300 K) 13 exp(— 7/T) and a temperature range at 
which the rate coefficients hold. For some reactions, more than 
one entry with non-overlapping temperature ranges and respec- 
tive fit parameters is present because of multiple sets of experi- 
mental data. In the given fit form, 7 is equivalent to the reaction's 
activation energy, the energy barrier that has to be overcome for 
the reaction to start. When fitting the measured reaction rates 
over a particular range of temperatures, a negative 7 may some- 
times give the best fit. Special care has to be taken when using 
these rate coefficients outside their given range. Especially at 
very low temperatures, which are relevant for many astrophysi- 
cal situations, negative 7 may lead to unphysically high reaction 
rates. 

To prevent unreasonable results, Woo dall et al.l (120071) give 
a list of reactions whose rate coefficients should be set to zero 
at 10 K. This may be feasible in chemical models, where the 
gas temperature is a constant, given parameter. However, in the 
framework of, e.g., numerical models of photon-dominated re- 
gions (PDRs), where the temperature is determined from de- 
tailed energy balance and may cover a wide range (~5 K- 
5000 K), it is preferable to use reaction rate coefficients, which 
do not sho w discontinuities across the applied temperatur e range 
(see. e.g.. lHoilenbach & TielenJl999URolhg et al.ll2007h . 

In Figure Q] we show two example reactions that have a di- 
verging reaction rate coefficient at low temperatures because of 
7 <C 0. Using these extrapolated values in a chemical network 
calculation may lead to artificial and very unrealistic results! 

To illustrate the problem, we compare the results of three 
PDR model calculations using the KOSMA-r PDR model 
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Table 1. All reactions with negative 7 coefficients that were refitted. The coefficients on the left are the original UDFA values, the 
new fits to the coefficients are given on the right side (numbers in parentheses mean: times ten raised to that power). The last column 
gives the relative error brackets across the given temperature range (given in K) with respect to the original fit. 
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Fig. 1. Two example reactions to demonstrate how extremely reactions might diverge from physically reasonable numbers if 7 <C 0. 
The final fit is shown in red, the UDFA rates are shown in light blue. Left The UDFA extrapolation to fc(10K) = 4.5 x 10 4 compared 
to 2.5 x 10~ 12 in our new fit. RightThe UDFA extrapolation to fc(10K) = 2.3 compared to 1.3 x 10~ 10 in our new fit. 



(Rollig et al. 2006). We calculated spherical model clouds with 
mass M = 1 M Q , surface gas density n = 10 6 cm~ 3 , and FUV 
illumination XDraine = 10 6 . The chemical network consists of 
198 species and 3237 reactions. Except for the surface forma- 
tion of molecular hydrogen, a pure gas-phase chemistry was cal- 
culated. The only difference between the models is the applied 
set of reaction rate coefficients. In the first model we used the 
original parametrization of the reaction rate coefficients^. In the 
second model we deactivated these reactions by setting their re- 
action rates to zero. In a third model we used our new fit to 
the reaction rate coefficients from Table [TJ In Figure |2]we show 
the relative sulfur abundance versus Ay for all three cases. For 
Ay > 7 the models start to deviate significantly. The relevant 



3 Using the original UDFA rates may lead to non-converging models 
depending on the given parameters. The model shown in Figure[2]is one 
of the few that converged at all within our parameter grid. 



reaction for this comparison is 

2 + S — >SO + . (1) 

The original UDFA rates are valid for T > 200 K (see Fig. [TJ 
left). The model temperature goes down to 40 K at Ay = 3.5. At 
this temperature, the UDFA reaction rate coefficient is a factor 
670 higher than at 200 K, leading to a very efficient destruc- 
tion of S. At Ay = 20 reaction Q] is three orders of magnitudes 
stronger than the next important destruction reaction. Switching 
off this reaction allows the atomic sulfur to remain in the gas up 
to large Ay ■ Our new rate coefficients give a result between these 
two extremes. In our comparison, the choice between the scenar- 
ios a), b), and c) strongly affects the whole sulfur chemistry at 
large Ay. 

The UDFA lists the accuracies of the given fits, usually rang- 
ing between 25-50%. However, outside their given temperature 
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Fig. 2. Relative sulfur abundance as function of Ay for a spheri- 
cal cloud. The model parameters are mass: M = 1 M Q , surface 
density n = 10 6 cm~ 3 , and FUV intensity x = 10 6 in units 
of Draine. Three different chemical networks were applied: a) 
the original UDFA rates, b) the original UDFA rates with all re- 
actions possessing 7 < set to zero, and c) all reactions with 
7 < replaced by the new fits from Table [TJ 



ranges, negative 7 may lead to errors many orders of magni- 
tude larger. For these reactions we calculated new fit parameters 
to allow for continuous rate coefficients even when T becomes 
very low. Furthermore, we calculated a single set of fit parame- 
ters for those reactions that have multiple entries in the database. 
These multiple sets may or may not have overlapping temper- 
ature ranges. Our final fit is designed to be valid over the total 
temperature interval with only small deviations. 

It is well possible that the results obtained with our new 
fits have large errors, especially outside the original tempera- 
ture regime. The peculiar shape of the fit in Figure fl]below T m j n 
is a direct result of the choice of the parametrization function 
in UDFA and it is unlikely that the low temperature behavior of 
this reaction follows that particular shape. However, the diver- 
gent low temperature rates of the original UDFA parametrization 
are unphysical, and the new fit is our best attempt at improving 
this behavior while minimizing any new uncertainties. 



2. Fit procedure 

Below we describe the applied fit strategy. Generally, we are 
willing to sacrifice some of the originally given accuracy to gain 
continuity down to very low temperatures. To be precise, we did 
not calculate a new fit to previously attained experimental data. 
We calculated a new fit to rate coefficients calculated from the 
original UDFA fits over the given range. The steps are the fol- 
lowing: 

1 . tabulate the reaction rate coefficients over the given temper- 
ature range, using the given fit-accuracy as error. These are 
our artificial data. 

2. if 7 <0, fixfc(T = 10K) 

3. calculate new fit parameters a, (3, ^ 

4. calculate the error of the new fit with respect to the original 
UDFA parametrization. 



The fitting was made with Mathematica from Wolfram Research, 
using a Levenberg-Marquardt algorithm. The fitting procedure is avail- 
able upon request. 



In the appendix we show plots of all new fits presented in this 
paper. 

2. 1. Discussion of the fit strategy 

The details of the fit strategy may have considerable influence on 
the final fit result. Particularly point 2 prevents the reaction rates 
from becoming unphysically high. Of course, this procedure in- 
troduces new, possibly very large errors to the rate coefficients. 
But these errors are usually of the same order as those from the 
UDFA fits. We now consider some major influences on the fit 
results. 

2.1 .1 . Placement of 'data' points 

The number and placement of 'data' points has a direct influ- 
ence on the quality of the fit because it implicitly introduces a 
weighting of the respective temperatures. Using the Levenberg- 
Marquardt algorithm is in effect applying a sequence of linear 
least-squares fits to a nonlinear function. Regions with more 
'data' points contribute more to the total sum of squares and 
force the fit function closer to these 'data' points. For instance, 
placing temperature points equidistant in log-space will place 
more grid points at lower temperatures compared to equidis- 
tant linear spacing, weighting lower temperatures stronger than 
higher temperatures. Optimal fit results could be achieved by 
adaptively changing the number and positions of grid points un- 
til the error is minimal, but to keep a consistent error estimate, 
we used the same temperature gridding for all fits. To fit reac- 
tions with 7 < 0, we employed an equidistant temperature grid 
in log-space with 20 grid points, while we used linear spacing 
for fitting reactions with multiple entries, using 20 grid points 
per temperature interval. 

2.1.2. Fixing k(T = 10K) 



C + O, — > CO + o 




Fig. 3. Comparison of two fit results for different choices of 
k(T= 10K). The UDFA rate coefficient, valid between 295 and 
8000 K, is shown in light blue. The best fit for k{10K) = 10^50 
is shown as red, dashed line (a — 1.01(— ll),/3 = 0.23,7 = 
913). The best fit for k(10K) = fc(T m ; n ) is shown as red, solid 
line (a = 9.95(-12), = 1.05, 7 = -42.5). 

Choosing a fixed k(T = 10K) prevents the reaction rates 
from becoming unphysically high. However, the exact choice of 
k(T — 10K) is completely arbitrary since no real measurement 
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is available. Moreover, it poses a major constraint to the fit and 
an unfortunate choice will le ad to very large devia tions from the 
true reaction rate coefficient. IWoodall et al.l C2007) suggest to set 
k(10K) = 0. It is possible that lim T _>.io k(T) = 0. It is also 
possible that k(10K) >> 0, we simply do not know. To mini- 
mize the deviations from the only available 'data', we chose to 
set fc(10A") = fc(T min ) if 7 < with T min the minimal valid 
temperature from the UDFA parametrization. This also reduces 
any over/under-swinging behavior of the new fit. To illustrate 
how the choice of k(10K) influences the form of the final fit 
function we compare in Figure [3]the fit for fc(10A") = fc(T m ; n ) 
(red, solid line) and for fc(10A) = 10~ 50 (red, dashed line). 
In this particular case, the latter leads to significantly stronger 
deviations from the UDFA parametrization than for k(10K ) = 

fc(Tmin)- 



mentfl Instead, whenever we are speaking of fit errors, we 
actually compare how two different models deviate from each 
other. Nevertheless, because our aim is to change the original 
parametrization as little as possible, we will use this deviation to 
quantify the quality of our fits. 

In Table Q] and [2] we give the relative error brack- 
ets of our new fit, i.e. [min(ei),max(ej)], with a = 
(kf(Ti) - k u (Ti)) /k u (T z ) for all T t e [T min , T max ]. Here, k u 
is the UDFA parametrization of the reaction, and k* is our new 
fit. 

2.2. Fitting reactions with negative 7. 



2.1.3. Calculating new fit parameters 




linear fit function 

logarithmic fit function 



6000 sooo 10000 

Temperature [K| 



Fig. 4. Comparison of fit results for a linear fit function k(T), 
shown as solid, red line, and a logarithmic fit function log(fc(T)), 
shown as dashed, red line. 



We fitted the fit function k(T) = a (^) e~T to the list of 
'data' points created in steps 1 and 2. In most cases it is better to 
perform the fit in log-space, i.e. use log(fc(T)) instead of k(T) as 
fit function. Otherwise, the Levenberg-Marquardt fit would have 
the tendency to fit k(T) to absolutely higher temperatures, while 
we are interested in the lower temperature regime. In some cases 
however, a linear fit function is more appropriate. This case is 
shown in Figure |4] We try to find a simultaneous fit across the 
multiple temperature regimes of the reaction C + N — > CN. The 
linear fit produces a continuously declining result for T — > 0, 
while a logarithmic fit function produces the opposite behavior, 
which results in large deviations of the 10 K rate predicted by the 
UMIST parametrization and our fit result. Again, the true reac- 
tion rate behavior at very low temperatures is unknown, hence, 
our aim is to minimize any deviation from the UDFA fit. The 
logarithmic fit differs from the 10 K UDFA rate by 8 orders of 
magnitude, while the linear fit differs only by a factor of 40. For 
this reaction we chose a linear fit function, for all other fits in 
this paper a logarithmic fit function was used. 

2.1.4. Error calculation 

Strictly speaking, the term error is misleading in our con- 
text, since it suggest model deviations from real measure- 




500 1000 1500 2000 

Temperature [K] 

Fig. 5. Example reaction that can be fitted very well (maximum 
additional error is 11%), but has a less extreme low temperature 
behavior. The final fit is shown in red, the UDFA rates for each 
individual temperature regime are shown in light blue. The dif- 
ference in k(10 K) is 10 orders of magnitudes between the UDFA 
rates and our new fitted rates. 



To obtain a new fit for reactions with 7 0, we used a log- 
spaced temperature grid with 20 grid points and a logarithmic fit 
function. Collider reactions, i.e. reactions of type "CL" in UDFA, 
were excluded from the fit. Apparently ist is possible to find a 
very good fit for some reactions without any divergence at low 
temperatures and without introducing any significant additional 
error. The prime example for this is shown in Figure Using the 
UDFA reaction rates gives an unreasonably high reaction rate 
coefficient at 10 K of 5! The new fitted rate coefficient at 10 K 
is 2 x 10~ 10 . Even though we do not know the real reaction 
strength at 10 K, we consider our new result to be much more 
reasonable compared to the extrapolated numbers from UDFA. 

Some of our new fits even have negative 7, which shows that 
7 < not always has to lead to divergent behavior. If 7 > —100, 
the rate inclination for low temperatures is moderate and should 
not lead to unreasonable results^. In Table Q] we present the fit 
results for reactions that are listed in UDFA with 7 < 0. We give 
the original and the new fit parameters as well as the relative 
error brackets of the new fit. 



The error bars in some of the figures show the uncertainties of the 
UDFA reaction rates. These errors are not taken into account by the 
Levenberg-Marquardt algorithm when fitting a new parametrization. 

8 UDFA list more reactions with 7 < than given in Table [T] We 
omitted those reactions because they showed no significant divergence 
at low temperatures. 
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Fig. 6. UDFA lists three C + recombination reactions for temper- 
atures from 10 to 41000 K. The two reaction rates for the higher 
temperatures show a very different behavior compared to the re- 
action that is valid from 10 to 8000 K. Applying different values 
of T max during the fit leads to significantly different results. 



He + + CO — > O + C + + He 
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Fig. 7. Example reaction with apparently mutually exclusive re- 
action rates listed in UDFA. The final fit is shown in red, the 
UDFA rates for each individual temperature regime are shown in 
light blue. The two temperature ranges overlap and the respec- 
tive fits to the reaction rates, here shown in light blue, display a 
very different behavior. 

Some of the reactions given in Table Q] also appear in 
Table [2] because they have counterparts with different tem- 
perature regimes. Some of the numerical problems described 
above can be circumvented by choosing an alternative reac- 
tion rate coefficient with a more appropriate temperature regime. 
Nevertheless, we list these reactions for cases where that is not 
possible. 

2.3. Fitting reactions with multiple entries. 

The fits to reactions with multiple database entries are, in most 
cases, much less sensitive to the details of the fit strategy, but, 
overall, the errors of the simultaneous fits are larger compared to 
the fits to reactions with 7 < 0. Here, the partly very inconsistent 
entries in UDFA for the different temperature ranges contribute 
the most to the errors in the final, simultaneous fit. For cases thta 
span a very large temperature range, the choice of T max when 
setting up the artificial 'data' points poses a major constraint on 
the fit result. This is demonstrated in Figure [6] 



H + NH, — > NH + H 2 



ur"< 




500 1000 1500 2000 2500 3000 

Temperature [K] 



Fig. 9. Simultaneously fitted reaction rates in this figure have a 
maximum relative error of 33%. The final fit is shown in red, the 
UDFA rates for each individual temperature regime are shown 
in light blue. 

In Table |2] we present the fit results for reactions with multi- 
ple entries in UDFA. The valid temperature ranges of these re- 
action rates may or may not overlap. Even worse, some of the 
multiple entries appear to be mutually exclusive. Especially the 
high temperature rate coefficients for the reaction He + +CO— 
+ C + + He appear to be at odds with both theory and experiment 
and it may be advisable to only use the low temperature rate co- 
efficients (see Figure|7]i. However, this assessment is far beyond 
the scope of this paper. We apply our fit routine to all reactions in 
UDFA with multiple entries and leave it to the informed reader 
to decide whether to use the simultaneous fit or whether to select 
a reaction with an appropriate temperature regime. 

For the fitting, we limited the upper temperature range to 
20000 K where T max is higher, thus neglecting higher tempera- 
tures. The errors are explicitly calculated for this range and may 
be significantly larger for higher temperatures. Despite the ap- 
parently large errors we consider the fits relatively good in most 
cases. Remember that we fit multiple, sometimes very disparate 
sets of 'data' points, and therefore, large errors are unavoidable. 
This is demonstrated in Fig.|8]for two cases from Table |2] 

The reader should also recall that the reaction rates from 
UDFA usually have errors on the order of 25-200%. One could 
be tempted not to use the UDFA rates to prevent the errors to rise 
even more because of our simultaneous fitting. However, by do- 
ing so one might implicitly accept much larger errors when mov- 
ing from one temperature regime to the other. This is demon- 
strated in Figure [9] 

The simultaneously fitted reaction rates in Figure |9] have a 
maximum relative error of 32%. This may appear large com- 
pared to the individual errors for each temperature regime. But 
when only using either one of the two reactions we would im- 
plicitly accept a very large error when the temperature is around 
400-500 K, where both individual rates disagree by one to two 
orders of magnitude. Our simultaneous fit is a reasonable fit in an 
attempt to keep the experimental findings for both temperature 
regimes. 

The special case of carbon recombination is shown in 
Figure [6] Three different sets of fit parameters are given in 
UDFA for three distinct temperature ranges up to 41000 K. 
Especially the high temperature rates show a very different be- 
havior compared to the low temperature range. A simultaneous 
fit gives very large errors because of the large deviations from the 
high temperature rates. However, for most applications, temper- 
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Fig. 8. Two example reactions to demonstrate the wide range of error values. The final fit is shown in red, the UDFA rates for each 
individual temperature regime are shown in light blue. Left This reaction is listed with two very differently behaving reaction rates. 
Consequently, the simultaneous fit has relatively large errors especially at the interface between bow temperature regimes. Still, 
we consider the fit to be very good because of its overall capability to reproduce the general behavior in both regimes. Right This 
reaction is listed with three distinct temperature regimes. The simultaneous fit is nevertheless a very good approximation to all three 
UDFA rates in the relevant temperature ranges. 



Table 2. All reactions with multiple reaction rate coefficients for several temperature ranges. The last column gives the relative error 
brackets. For the fitting, we limited the upper temperature range to 20000 K, neglecting higher temperatures. The errors for higher 
temperatures might be significantly larger. The temperature limits for UDFA are given in Kelvin. 





Educts 






Products 






Tmin 


Tmax 




/^ncw 




[error] 
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+ 


NH 2 


— > NH 


+ 


H 2 






73 


3000 


4.61(-12) 


1.02 


2161 


-0.05, 0.32 


C 


+ 


o 2 


— > CO 


+ 


O 






10 


8000 


5.56(-ll) 


0.41 


-26.87 


-0.50, 2.89 


CH 


+ 


O 


— ► CO 


+ 


H 






10 


6000 


6.02(-ll) 


0.10 


-4.46 


-0.17, 0.13 


N 


+ 


OH 


— ► NO 


+ 


H 






10 


2500 


6.65(-ll) 


-0.23 


0.24 


-0.12, 0.26 


O 


+ 


OH 


-> o 2 


+ 


H 






10 
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5.28(-ll) 


-0.46 


4.14 


-0.22, 0.66 


O 


+ 


c 2 


— ► CO 


+ 


C 






10 


8000 


1.25(-10) 


0.47 


-16.13 


-0.27, 0.32 


H 2 


+ 


CH 4 + - 


-> CH 5 + 


+ 


H 






15 


41000 


4.89(-ll) 


-0.14 


-36.14 


-0.47, 0.96 


H 2 


+ 


NH 3 + 


— ► NH 4 + 


+ 


H 






10 


41000 


3.09(-13) 


1.08 


-50.87 


-0.80, 0.83 
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+ 


CO 


— ► O 


+ 


C+ 


+ 


He 


10 
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2.27(-10) 
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-0.83, 1.51 


CH 3 + 


+ 
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+ 
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10 
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-0.40, 0.36 


o+ 
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+ 
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23 


41000 


2.42(-12) 
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-0.19, 1.61 
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+ 
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+ 


H+ 






10 
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8.41(-16) 


0.43 


-9.59 


-0.36, 0.30 


H+ 


+ 


O 


-> o+ 


+ 


H 






10 


41000 


6.86(-10) 


0.26 


224 


-0.08, 0.05 


NH 2 + 


+ 


e~ 


— ► N 


+ 


H 


+ 


H 


12 


12400 


1.78(-7) 


-0.80 


17.14 


-0.23,0.31 


NH 2 + 


+ 


e _ 


— i NH 


+ 


H 






12 


12400 


9.21(-8) 


-0.79 


17.11 


-0.23,0.31 


NH 4 + 


+ 


e — 


-+ NH 2 


+ 


H 


+ 


H 


12 


34670 


4.67(-7) 


-1.25 


41.88 


-0.50, 4.02 


NH 4 + 


+ 


e~ 


NH 2 


+ 
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-1.25 


41.91 


-0.50, 4.02 


NH 4 + 


+ 


e — 
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+ 


H 
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1.54(-6) 


-1.25 


41.92 


-0.50, 4.02 
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e~ 


—> H- 


+ 
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3.37(-16) 


0.64 


9.17 
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C+ 


+ 


e — 


-> C 


+ 


hv 






10 
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2.36(-12) 


-0.29 


-17.55 


-0.58, 0.41 


N+ 


+ 


e~ 


— ► N 


+ 
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10 


41000 


3.5(-12) 


-0.53 


-3.18 


-0.34, 0.08 


H+ 


+ 


H 


-> H 2 + 


+ 


hu 
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32000 


1.15(-18) 


1.49 


228 


-0.18, 0.42 


C 


+ 


N 


— ► CN 


+ 


hv 






300 


14700 


5.72(-19) 


0.37 


50.95 


-0.16, 0.46 


C 


+ 


O 


-* CO 


+ 


hv 






10 


14700 


4.69(-19) 


1.52 


-50.51 


-0.69,3.10 


c+ 


+ 


O 


— ► CO+ 


+ 


hv 






10 


13900 


2.44(-18) 


-0.07 


2.77 


-0.09, 0.06 


o 


+ 


Si+ 


— ► SiO+ 


+ 


hv 






10 


14700 


5.92(-18) 


-0.67 


19.37 


-0.45, 4.87 


H 


+ 


HNC 


— > HCN 


+ 


H 






10 


1000 


1.14(-13) 


4.23 


-115 


-0.70,5.10 



atures up to 8000 K are sufficient and we can neglect the errors 
for high temperatures. In that case, the new fit has much smaller 
errors. For many practical purposes, one could even completely 
discard the two higher temperature ranges and only use the low 
temperature rates. For the sake of completeness, we still list the 
new fit temperatures up to 20000 K. 

2.4. Discussion 

For three reactions in Table [T]and|2]one finds updated measure- 
ments and theoretical values of the reaction rate coefficients in 



KIDA k K (T). A detailed comparison of the respective reaction 
rate coefficient behavior is beyond the scope of this paper. We 
will only take a look at the updated recommended 10 K val- 
ues for these reactions and compare them to our new values 
^(lO^inTableEJ 

For all three reactions, the 10 K rate coefficients achieved 
from our new fits are within a factor of two compared to 
the recommendations from the KIDA Kinetic Database for 
Astrochemistry. For the reaction O + OH, our value is only 40% 
higher than the recommended value. Given the overall uncertain- 
ties, the agreement between our new fits and the KIDA values is 
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Table 3. Comparison of the recommended 10 K reaction rate coefficients available in KIDA k (10K) with k' (10K) from our new 
fits. 





Educts 


Products 




k K (WK) 

^ — i 
cirrs 


kl(lQK) 

°. — 1 
cirrs 




References 


N 


+ OH - 


-> NO + 


H 


4.5 x 1CT 11 


8.44 X 10" 


-11 


Jorfi et al. (2009) 


O 


+ OH 


-> H + 


Oz 


4 x 10" 11 


5.57 x 1CT 


-11 


Liaue et al. (2009); Ouemener et al. (2009) 





+ NH 2 - 


-> HNO + 


H 


1 x 10" 10 


4.80 x 10" 


-11 


Baulch et al. (2005) 



remarkable and an indication that our general fit approach is vi- 
able to achieve reasonable results. 



3. Summary and conclusions 

We presented an approach to overcome some numerical and 
practical difficulties in present-day databases of chemical reac- 
tion rates. One of these difficulties is the existence of fits to re- 
action rate coefficient that unphysically diverge at very low tem- 
peratures. This is especially problematic for numerical codes, 
which calculate chemical networks down to very low tempera- 
tures and thus have to use many of the reaction rate coefficients 
that are listed in chemical databases beyond their given temper- 
ature regime. We demonstrate that it is possible to find a new 
parametrization of these reactions that ensures a less erroneous 
behavior and allows calculations to temperatures of 10 K and 
lower. 

A second complication that occurs in chemical databases is 
that some of the listed reactions have multiple entries because 
the experimental data were obtained for different temperature 
ranges. Any user is instructed to use the listed reaction rate co- 
efficients only in their valid temperature regime. Unfortunately, 
multiple entries in these databases very often display fairly in- 
consistent properties. This is particularly problematic for numer- 
ical codes that continuously change the temperatures when cal- 
culating chemical networks because it introduces artificial nu- 
merical jumps in the coefficient matrix of the equation system to 
be solved. One solution could be to restrict oneself to only use 
the one reaction rate that has a temperature regime closest to the 
computational regime. However, in many cases this is not pos- 
sible because the parameter space relevant for application spans 
multiple temperature regimes. Additionally, this means that one 
implicitely accepts quite large errors when leaving the valid tem- 
perature regime. We suggest to replace these multiple entries 
with a single, new rate which is fitted simultaneously to the in- 
dividual ones. This minimizes the error when moving from one 
temperature regime to the other and additionally overcomes the 
problem of discontinuous reaction rate jumps. 

We consider the new parametrization presented in this paper 
only a temporary solution until better experimental data for low 
temperature reaction rates are available. 

As a side effect we demonstrate that in some cases one might 
find an equally good parametrization of the reaction rate coeffi- 
cients without any artificial numerical side effects just by taking 
care not to allow too low negative values of 7. The chemical 
databases try to express the possibly very complicated physics 
of the various chemical reactions with only three fitting param- 
eters. Of course, this is an extremely strong constraint on any 
parameter fit to the complex experimental data. Although, un- 
der certain circumstances, chemical reactions listed in chemical 
databases like UDFA, OSU, and others might reasonably be ap- 
plied even outside their given scope, it is important to ensure that 
numerical side effects from parametrization are minimized when 



possible. For many numerical applications this is much more im- 
portant than a smaller total error of the reaction rate coefficients 
in a limited range in temperature. 

Any user should always take care not to blindly use reaction 
rate coefficients from chemical databases outside their tempera- 
ture scope. This is also true for the new parametrizations listed 
in this paper. In each case, one has to balance advantages against 
disadvantages. If you are sure that you will keep within the given 
temperature regime, then there is no need to use our rate coeffi- 
cients because they add an additional error to the computation. 
Under certain circumstances, this additional error might be neg- 
ligible compared to the numerical problems that come with the 
original rate coefficients. 
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Appendix A: Plots of all new fits. 

Here, we show the plots of all new fits to allow the user an assess- 
ment of the fit quality over the entire fitted temperature range. 
The plots are shown in the same order as in Table [T]and|2] 
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C + NH 2 — * HCN + H C + NH : — » HNC + H C + 2 — > CO + O 
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Fig.A.l. 



Fits to reactions with negative 7. 




Fig. A.2. Fits to reactions with negative 7 (continued). 
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Fig. A.3. Fits to reactions with multiple entries. 
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Fig. A.4. Fits to reactions with multiple entries (continued). 



